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ABSTRACT 

Helium star-compact object binaries, and helium star-neutron star binaries in par- 
ticular, are widely believed to be the progenitors of the observed double neutron star 
systems. In these, the second neutron star is presumed to be the compact remnant of 
the helium star supernova. Here, the observational implications of such a supernova 
are discussed, and in particular are explored as a candidate 7-ray burst mechanism. 
In this scenario the supernova results in a transient period of rapid accretion onto 
the compact object, extracting via magnetic torques its rotational energy at highly 
super-Eddington luminosities in the form of a narrowly beamed, strongly electromag- 
netically dominated jet. Compton scattering of supernova photons advected within 
the ejecta, and photons originating at shocks driven into the ejecta by the jet, will 
cool the jet and can produce the observed prompt emission characteristics, including 
the peak-inferred isotropic energy relation, X-ray flash characteristics, subpulse light 
curves, energy dependent time lags and subpulse broadening, and late time spectral 
softening. The duration of the burst is limited by the rate of Compton cooling of the 
jet, eventually creating an optically thick, moderately relativistically expanding fire- 
ball which can produce the afterglow emission. If the black hole or neutron star stays 
bound to a compact remnant, late term light curve variability may be observed as in 
SN 2003dh. 
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1 INTRODUCTION 

Despite beinR d iscovered more than 30 years ago 
iKlebesadel et aDirgySjl , and the fact that bursts have been 
detected daily since the advent of the Compton Gamma-Ray 
Observatory in 1991, there still does not exist a satisfactory 
theoretical model which encompasses the entire 7-ray burst 
event. Nonetheless, in spite of their strong heterogeneity, the 
abundance of observations has led to a well described burst 
phenomenology, and is briefly summarised below. 

Burst durations are bimodal, with long bursts last ing 
~ 10 s while short bursts last ~ Is (Kouvcliotou ot all 
Il993^ . Long bursts typically have inferred isotropic luminosi- 



' erg/s (s 



itro pic lummosi- 
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ties on the order of ~ 10 
and references therein), although achromatic breaks in the 
afterglow light curves suggest strong colhmation (with open- 
mg ang les ~ 5-10°. iFrail et al. 2001) which, in turn, implies 
actual luminosities on the order of ^ io**~^^ erg/s. The ma- 
jority of the prompt emission is in the form of apparently 
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nonthermal 7-rays, well fit by a broken power law, typically 
peaking near 100-1000 keV, and softening throughout the 
burst. The prompt emission is composed of a large number 
of subpulses with typical widths on the order of a second in 
which lower energy emission lag behind, and are wider than, 
the higher energy emission. 

Long bursts are followed by optical and radio after- 
glows thought to be generated by a hot fireball interact- 
ing with the interste llar medium (see, e.g.. iMeszarosI 120021 : 
iLi fc Chevaiieill2003D . These can last many weeks and have 
in at least two cases (GRB 980425 and GRB 030329) in- 
cluded type Ib,c supernova (SN 1998bw and SN 2003dh, re- 
spectively) light curves 7-10 days after the prompt emission. 
The total bolometric energy in these bursts has been esti- 
mated from late time radio observ ations to be ~ 10^ ^ erg 
and roughly constant among bursts iBerger et al ] l2004ll . To 
date, there have been no instances of afterglow or supernova 
being associated with short bursts. 

Due to similarities in their temporal structure, dura- 
tion, and spectra. X-ray flashes appear to be low luminosity 
cousins of 7-ray bursts. Recently it has been shown that their 
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bolometric luminosity is indeed comparable to normal bursts 
jSoderberg et al .12004) . Despite this, their inferred isotropic 
energy is substantially less 10^^ erg) placing X-ray flashes 
up on the peak-inferred isotropic energy relation discovered 
by lAmati et al.l ll2002fl . and supporting a unified interpre- 
tatio n of X-ray flashes and 7-ray bursts ijSakamoto et alJ 
I2OO4I) . 

A number of mechanisms have been suggested for 7- 
ray bursts. Frequently, these address either the central en- 
gine and emission separately. Central engine models in- 
clude black hole and/or neutron star collisions (Eichlcr ct al. 
ll989l:IPac zvnski 1991; Nar avan et al.lll99i). maeneta r birth 
JUsovlll992l:iThompsoniil994l:lThompson et aljl2004l) . black 
hole birth JVietri fc Stellalll998^ ■ and coUapsar (see, e.g., 
iMacFadven et al.l 120011)" models. In most of these it is not 
clear how the spectral and temporal structure of the burst 
is produced. In contrast, there are also a number of mod- 
els which focus upon the e mission. These include the can- 
nonball l|Dado et al\ |2 002|. a nd references the rein), shot 
gun f Heinz fc BeEelman'l999), Comptonised jet (Thompson 
11994:. Ghisellini et al. 2 000'). and internal shock ( Piran 2004. ) 
models. Despite their ability to reproduce many of the ob- 
served spectral and temporal features of the bursts, few of 
these address directly the power source of the burst. As a re- 
sult, currently there are few models which can explain both 
the properties of the observed emission and the prodigious 
energy of observed bursts in a satisfactory manner. 

Here a model involving supernovae in helium star-black 
hole and helium star-neutron star binaries is presented. In 
this scenario the supernova results in a transient period of 
rapid accretion onto the compact object, extracting via mag- 
netic torques its rotational energy at highly super-Eddington 
luminosities in the form of a narrowly beamed, strongly elec- 
tromagnetically dominated jet. The prompt emission is pro- 
duced by Compton scattering supernova photons advected 
within the ejecta and photons created at shocks driven into 
the ejecta by the jet. The duration of the burst is limited by 
the rate of Compton cooling of the jet, eventually creating an 
optically thick, moderately relativistically expanding fireball 
which can produce the afterglow emission. In comparison 
with the previous models this has two advantages: firstly 
it provides a unified model for the central engine and the 
prompt emission, and secondly, since helium star-neutron 
star binaries are widely believed to be the progenitors of 
the observed double pulsars, it is assured that supernova in 
these systems will occur in sufficient quantity (see section 
13.51 . For convenience, typical values for the pertinent quan- 
tities which describe the model are collected in Tableland 
a schematic of the process is shown in Figure Q 

The formation and energetics of the jet are discussed 
in Section |5| Sections |3 and 2] describe the mechanisms by 
which the kinetic energy of the jet is converted in the the 
observed prompt emission and subsequent afterglow, includ- 
ing the limits this places upon the jet dynamics. Expected 
observational implications of this model are discussed in Sec- 
tion |K| Finally, concluding remarks are contained in Section 

m 
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pan 




Supernova Ejecta Velocity 


Van 


10^ km/s 


Neutron Star Mass 
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1.4 M0 


Neutron Star Radius 


Rns 


10'5 cm 


Black Hole Mass 




10 Mq 


Jet Lorentz Factor 


r 


20 


Jet Plasma a 


cr 


10^ 



Table 1. A number of model parameters are listed together with 
some of the jet characteristics obtained in the following sections. 



2 JET FORMATION 

A high power jet (or highly coUimated relativistic outflow) is 
a central feature of most 7-ray burst models. Observational 
constraints require that the luminosity of a jet with opening 
angle ~ F"^ (where F is the bulk Lorentz factor) be on the 
order of 

^--4?^^--^x^0^^(i)"(l(fe) -^/^' 

(1) 

where Li^o is the isotropic equivalent luminosity of the burst. 
For a burst of duration T, this corresponds to a total energy 
budget of 

6...6xlO"(i)"(j5i^)(^)»,.. W 

The observational features of the model presented here 
depend only upon the existence of a suitably luminous and 
electromagnetically pure jet (see, e.g., section f3.2t in the 
vicinity of the supernova, and are otherwise independent of 
the mechanism by which such a jet is produced. Nonethe- 
less, for concretness, here the jet is presumed to be the result 
of accretion driven electromagnetic torques upon a central 
compact object. Therefore, possible mechanisms for the pro- 
duction of such a jet are discussed below in the cases where 
the compact object is a black hole or a neutron star. How- 
ever, this is in no means meant to be an exhaustive discus- 
sion of the problem of jet formation, which is beyond the 
scope of this paper. 

2.1 Black Hole Binaries 

A considerable literature exists regarding the produc- 
tion of jets via the electromagnetic extraction of an- 
gular momentum from ac c reting black holes (se e , e.g. , 
Blandford & Znaiek 1977; 'Mc Kinnev fc Gammij l2004l : 
Komissarov. 2004: .LevinsQn. .2QQh . and references therein). 
The maximum energy that may be extracted by this method 
is determined by the irreducible mass, Mir: 

Etot = Mbhc' (1 - 

\ Mbh J 

1 .2 , J- 2 

-lO^Yerg. (3) 
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Figure 1. A cartoon of the mechanism is shown (not to scale) with the principle features, supernova ejecta, accretion driven jets, 
up-scattered entrained and shock produced photons from the ejecta, and the regions where the jet is axially optically thin (r ~ 10^^ cm) 
and becoming nonrelativistic. 



where j is the dimensionless angular momentum of the black 
hole, and the expansion is appropriate for j <^ 1 (see, e.g., 
iMisner et al. 1973). Comparing this to the energy required 
to power a burst (cf. equation |5J requires j > 0.01, and thus 
puts a rather weak limit upon the spin of the black hole. 

In the model presented here, the torque upon the black 
hole is due to a substantially super-Eddington accretion flow 
resulting from the accretion of the supernova ejecta. An es- 
timate for the mass accretion rate in this scenario is simply 
given by the Bondi-Hoyle rate. 



A> 2 A f GMbh 

M ~ 7rr-g^/9sn«sn — 47rpsn I ^ 

~ 2 X 10'^ g/s , 



(4) 



where the Bondi-Hoyle radius is given by 

2GMbh 



ran 



^sn ~l~ 



(5) 



in which the sound speed (cs) may typically be ignored com- 
pared to the ejecta bulk velocity. The resulting density in the 
vicinity of the horizon is then 



47rr^c V c 



Psn : 



(6) 



where is the horizon radius. This implies a magnetic field 
strength on the order of 



~ 2 X lO" G , 



(7) 
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presum ably created via the m aKneto-rotational instability 
(MRI) fHawley & Balbus"l99j). The rotational energy may 
then be electror naKnetically tapped by, e.g ., the Blandford- 
Znajek process jBlandford fc Znaieklfl977f) . Recent numer- 
ical simulations of accreting black holes have produced 
Poynting flux dominated jets with typical luminosities of 



f^RnSRH — 



1 /GMbh\ 



2x lO^Verg/s, 



(8) 



where SIbh = jc/2rh is a measu re of the angular velocity of 
the black hole jMcKinnevl2005l) . Therefore, rapidly rotating 
black hol es (?' > 0. 8 for the black hole parameters given here, 
cf. Gam mie et alJ i2004'l are easily capable of producing a 
Poynting dominated jet of sufficient luminosity. 

This scenario is similar to t he central engine m odels 
for the coUapsar scheme (see, e.g., |PoBham_et_ al.ll ll998ll . the 
primary difference being the accretion rate. In both cases, 
the accretion flow is neutrino cooled, and thus able to reach 
substantially super-Eddington (cf. with 10^* g/s) accretion 
rates. However, whereas the typical coUapsar accretion rates 
are on the order of IMq s~^, here they are a considerably 
lower IQ-^Mos-^ 



2.2 Neutron Star Binaries 

A second, and perhaps more speculative, mechanism by 
which the jet might be formed involves a rapidly rotating 
neutron star. As with the black hole, in this scenario the en- 
ergy is provided by the spin of the star. The total rotational 
energy available is 



Eto 



lO^^Lj^rg, 



-Mns-Rns^ns 
5 



-'NS 2 
-UJ 



5J?NS 



(9) 



where lu is the angular velocity in units of the breakup ve- 
locity. In order to be sufficient to power a burst this already 
requires that uj > 0.03 and hence the neutron star must be 
rapidly rotating (P < 15 ms) consistent with observations of 
millisecond pulsars. 

As with the black hole, the accretion of the supernova 
ejecta may mediate the extraction of the rotational energy. 
In the neutron star case, the Bondi-Hoyle accretion rate is 
given approximately by 



M ~ 5 X 10^^ g/s , 



(10) 



at whi ch, as i n the black h ole case, neu trino cooling domi- 
nates llChevalier 1989.: Fryer et al]|l996ll . Near the neutron 
star, this results in a density on the order of 



M 



10 Pr gem , 



(11) 



where Vr = PrC is the infall velocity near the neutron star 
surface. This implies an accretion magnetic field strength on 
the order of 



Bns 4 X lO'^A 



12^-1/2, 



(12) 



Note that this is a lower limit as the radial infall velocity 
may be substantially less than the speed of light. This field 
will reconnect with the native stell ar field on time scale s 
comparable to the infall time scale jGhosh fc Lamblll979tl . 



Therefore, it may not be necessary for the native field of the 
neutron star to initially be dynamically significant for it to 
subsequently strongly couple the star to the disk. 

The resulting spin-down torque has been studied pri- 
mar ily in the context of propellor a ccretion flo ws (see, 
e.g..lGhosh fc Lambll979l : lEksi et al 12005: Romano va et alJ 
l2004h . Unfortunately, the magnitude of the spin down- 
torque is not a settled issue. Nonetheless, recent numerical 
and analytical efforts have implied that 



M ~ bIrW 



(13) 



where Rd is the radius at which the magnetic field is domi- 
nated by the inertia of the disk, is the magnetic field at 
this radius, Wesc is the escape velocity from this radius, and 
r], the so-called fcistness parameter, is the rotation frequency 
of the neutron star in units of the Kep lerian velocity at the 
dissipation radius faomanov a et aLll20 04: Eksi ct al. 2005i). 
Clearly the minimum dissipation radius is the radius of the 
neutron star, and thus the maximum possible luminosity is 



-J/NSmax — Bj^<jRj^<jCLO (1 



(14) 



which is sufficient to drive a burst provided that B > 10^^ G 
and LO is near unity (P ~ 1ms). The former is satisfied if 
/3r — 0.1. The latter must be satisfied for a propellor to op- 
erate if the dissipation radius is near the surface, placing a 
more stringent constraint upon the spin of the star. It should 
be noted that such a high spin may be expected in neutron 
star-helium star binaries if the helium star's hydrogen en- 
velope was expelled due to accretion by the neutron star or 
common envelope evoluti on, as is commonly th ought to be 
the case in such systems jBethe fc Brownlll99^ . 

In the standard propeller theory the dissipation radius 
is assumed to be approximately the Alfven radius, where 
the pressure from the neutron star's native field balances 
the accretion ram pressure. However, in the hyperaccret- 
ing millisecond pulsar case (considered here) the dominant 
magnetic field is the accretion fiow field. Therefore, this case 
may be more analogous to accreting black holes (in which 
the dissipation radius is roughly that of the horizon) than 
to accreting pulsars. 

The manner in which the dissipated energy is trans- 
ported out of the system is not clear a priori. Numerical 
simulations of propellers have found the development of a 
magnetic chimney, suggestive of the development of a jet, 
suggesting that so me fraction of the energy leaves in a Poynt- 
ing fiux ( Romanov^^L||200^. This interpretation is sup- 
ported by the obvious parallels with the accreting black hole 
case, in which a significant fraction of the spin-down energy 
is channeled into a Poynting jet. However, a definitive an- 
swer will depend upon the details of the accretion fiow and 
the magnetic coupling to the star, and will likely have to 
await detailed numerical simulations. 



2.3 Jet Emergence 

In order to have observational consequences it is necessary 
for the jet to emerge from the growing supernova. This is 
analogous to th e problem of jet emerge nce in the coUapsar 
model (see, e.g.. lMacFadven et all 200l) with a reduction of 
the external medium density by eight orders of magnitude. 
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The ram pressure of the jet, 
'"^ Syrr^c 12-Kr'^c 

^15 / Lis. 



(15) 



1 X 10' 



/ 3 

erg/ cm 



105ierg/sy V 1012 cm/ 

can be compared to the typical pressures in the supernova 
ejecta, 

Psn ^ IpsnvL =i 2 X lO" erg/cm^ , (16) 
6 

(the thermal photon pressure will be comparable at the or- 
bital separations considered here) and hence the jet will eas- 
ily escape the ejecta. 



3 PROMPT EMISSION 

The presence of a jet alone is insufficient to produce a 7-ray 
burst. Also required is a mechanism by which the consid- 
erable kinetic energy flux of the jet can be converted into 
the observed prompt emission. Many such mechanisms have 
been discussed in the literature. However, for highly beamed 
emission, a minimum requirement is that the jet be opti- 
cally thin along the jet axis. This may be accomplished in 
a number of ways, including clumpy jets. However, in the 
context of the model considered here, in which a large num- 
ber of seed photons are available entrained in the supernova 
ejecta, this immediately suggests inverse-Compton scatter- 
ing as the prompt emission mechanism. This has the con- 
siderable advantage over internal shock scenarios of being 
capable of converting the kinetic energy of the jet into the 
prompt emission at efficiencies approaching unity (see, e.g., 
iLazzati et aI]ll999^ . 



3.1 Photon Collimation 

Due to the high bulk Lorentz factor, the scattered seed pho- 
tons will natu rally be coUimated to wi thin r"^ of the jet 
axis (see, e.g.. iBegelman fc Sikoralll987f) . However, in addi- 
tion, there are optical depth effects which will also serve to 
beam the scattered photons. 

In the rest frame of the jet electrons, the Thomson 
depth is given by 



T — GTUeCAt' . 



When transformed into the lab frame this gives 
r= I (jTUecV [1 - P cos, 6) At 



(17) 



(18) 



= j (jTUeT {I- P cos e)dr , 

where in both cases rie is the proper electron number density. 
For the two limiting cases of across the jet [6 ^ F"^) and 
along the jet {6 < F~^), the optical depth is 



2r/F was used. Therefore, despite the considerable difference 
in scale length, the optical depth along the jet is a factor of 
F less than that across the jet. As a result, it is possible to 
have a jet which is optically thin to photons within ~ F~^ 
of the jet axis while being optically thick to all others. This 
provides an additional coUimating mechanism and explains 
why it is possible to Compton scatter a large number of the 
seed photons incident on the jet while allowing the scattered 
burst photons to escape. 



3.2 Implications for Jet Type 

The requirement that at interesting radii ry < 1 limits the 
baryon loading of the jet. The power associated with ions in 
the jet is 



-L/ions — -p TleTTlpC 



10-* i ^ 



2iTr 
7 

GT 



(20) 



20 



10" cm 



Hence, the condition that the jet be optically thin along it's 
axis at r ~ lO^'^ cm requires that in the jet. 



- 1 > 10* , 



(21) 



where a is the ratio the Poynting flux to the kinetic energy 
flux of the baryons. Therefore, the jet must be extremely 
electromagnetically dominated. It should be noted however 
that observational precedent for such an electromagnetically 
dominated outflow exists in the context of the Crab pulsar 
in which cr 10^ a t these radii (see, e.g., IVlahaki j l2004l : 
iKomissarov fc Lvubarsk v 20o3) . 

The magnetic fields within the jet required to generate a 
Poynting fiux with the 7-ray burst luminosities will typically 
be of order 



B '. 



, Lisa 



3 X 10^ 



Lis, 



10^-*^ erg/s 



1/2 



10" cm 



(22) 



for r ^ -Rns^- The resulting synchrotron cooling time for 
an electron with Lorentz factor 7 in the frame of the jet is 
then 



ts 



10 



DTT ;i^7 



syn 
10 



cB^ 



(23) 



10^1 erg/s 



10" cm 



hence these electrons may always be treated as cold in this 
frame. 



ry = —GTUer 

Tj_ = 2aTner 



= 



TV 

2 ' 



(19) 



where Ue (x r ^ (changing the power law index changes ry 
by factors of order unity) was assumed and a jet width of 



1 The conical structure of the jet will not necessarily extrapolate 
down to the acceleration region which is expected to occur on a 
scale of many gravitational radii. If it were this would suggest 
magnetic field strengths on the order of 10^^ G, as is suggested 
in magnetar models of 7-ray bursts as opposed to the far more 
conservative few x 10^^ G discussed here. 
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3.3 Compton Seed Photons 

From momentum conservation, the energy of a photon after 
a single Compton scatter is given by 

1 — /3 cos 9i 



1 -/3 cos 61/ + e, (1 - cose) /Vm^c^ 



(24) 



where 9ij are the initial and final photon propagation an- 
gle with respect to the jet axis, and O is the angle between 
the initial and final photon propagation direction. For pho- 
tons scattered to within F"^ (as expected from relativistic 
beaming and optical depth coUimation) , the resulting energy 
is given approximately by 



£/ = min (2F^ei, FrrieC^) 



(25) 



Therefore, in order to generate a burst with an 
isotropic luminosity of Liso, a luminosity of Liao/SF* ~ 
10''^(F/20)~* erg/s of seed photons are required to impinge 
upon the jet. 

There are two sources of seed photons, the thermal pho- 
tons from the supernova itself and those produced by strong 
shocks driven into the ejecta by the jet. In the first case, due 
to the ejecta's high Thomson depth, these will be entrained 
in, and in thermal equilibrium with, the ejecta. Thus, assum- 
ing the ejecta cools adiabatically as it expands, the photon 
temperature will be 



3fc 



r 



keV, 



(26) 



and hence the typical photon energy is on the order of 1 keV 
for an orbital separation of > 2i?Hc which would be expected 
if the Helium star is filling its Roche lobe. The associated 
energy density of photons is 



ut = aT^ 



2 X 10' 



r 



lO" erg/cm^ , (27) 



(where here a is the Stefan-Boltzmann constant and not 
to be confused with the orbital separation) and thus the 
luminosity of thermal seed photons available to the jet is 
approximately 



2r' 



10* 



10^2 cm 



erg/s, (28) 



where clearly this depends upon the alignment of the jet 
relative to the orbital plane through the dependence of ut 
upon the jet position. In the second case, the kinetic energy 
of the ejecta is thermalised at the jet boundary providing a 
seed photon luminosity of roughly 



Ls 



1 3 2rf 

2 Psnf sn p 



10" 



(29) 



1012 cm 



erg/s , 



which is comparable to Lt- Depending upon the shock struc- 
ture and location these may also be thermalised. Due to its 



weaker dependence upon r {pa 



) this contribution 



will dominate at large radii. The total luminosity of seed 
photons entering the jet is then simply 



— Lt + Ls ■ 



(30) 



which is sufficient. 



3.4 Implications for Jet Lorentz Factor 

Recent simulations suggest that in vacuum it would be ex- 
pected that the jet would accelerate due to magnetic stresses 
until it was no longer electro magnetically dominated, i. e., 
F would approach a ~ 10'' iSpitkovskv fc Aronj |2004) ^. 
If the seed photons are indeed from a thermal distribution 
with T-y ~ 1 keV as argued in the previous sections, the up- 
scattered emission would peak at a few GeV, far higher than 
is observed. However, in the situation under consideration 
here, the Compton scattering will induce a drag on the jet, 
substantially reducing its velocity. This will manifest itself 
differently depending upon whether or not the jet is opti- 
cally thin along its axis. 

When it is optically thin, the rate at which energy is 
lost by the jet to the up-scattered photons at a given radius 
is simply given by the rate at which seed photons enter the 
jet and their subsequent up-scattered energy, i.e.. 



d-L _ disced 

dr dr ' 

where Lsecd as a function of radius is defined by 

2r 



-'^sccd('") = / {ut + us) Vsn- 
J rn 



-dr , 



(31) 



(32) 



and us is the energy density of seed photons produced in the 
shocks (of order psnWsn/2). Since dL/dr is negative definite 
it will generally produce a deceleration in the jet. 

In contrast, when ry > 1 the up-scattered photons can- 
not escape. They will consequently rescatter, introducing a 
second term to the energy loss: 



_ _ 2d£;secd „, dl^ 
dr dr dr 



= -A2F^i_,. 
dr 



(33) 



Unlike the optically thin case, for strongly decelerating jets 
this can produce a positive accelerating force. This provides 
an efi'ective mass loading of the jet, and since a is so high, 
will determine the maximum F reached, i.e., 



ijc 



-^iso 



2F2Ls. 



8F^ (ut + us) VanT^ 



30 



\ 10^-*^ erg/s 



1/4 



10" cm 



-1/2 



(34) 



Hence if the jet becomes optically thin near r = 2 x 10^^ cm, 
F ~ 20 as assumed thus far. After ry drops below unity 
the jet will continue to decelerate until it exits the ejecta or 
stalls. 

The Lorentz factors considered here are substantially 
lower than those implied by the typical resolution to the 
"compactness problem " (which are > few x 10^, see, e.g., 
iLithwick fc Sarill200 j) . The runaway pair production asso- 
ciated with the "compactness problem" is not present here 



2 Note that this is in contrast to many of the theoretical 
estimates in the literature in which behaviour expecte d by 



see. 



Le^ I199-: 



■ expecteg Dy 
lMicheJll969t 
; iBeskin et al.l 



one-dimensional models i s F — > a^ / ^ 
[Goldreich fc Ju UanI 1197(1 iBcgclman 

[l998). Spitkovs kv &: Arons 12004) claim that the discrepancy is 
largely due to the various artificial assumptions that these early 
efforts necessarily adopted, including spherical symmetry, split- 
monopole field configurations, and the region of applicability of 
the force-free approximation. 
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as a result of the low energy of the seed photons them- 
selves. In the jet frame, the seed photons have energy ~ Te^n, 
where ean is the typical seed photon energy in the lab frame. 
Since in the jet frame the electrons are cold, and thus all 
of the photon scattering is elastic, the maximum energy in 
any photon-photon collision is 2resn, which for F ~ 20 and 
Esn ~ 1 keV is insufficient to pair produce. 

Nonetheless, for a thermal seed photon distribution a 
high energy tail will exist, including photons with e > 
rrieC^/r and hence will be able to pair produce. The number 
density of such photons accumulated in the jet frame by the 
end of the burst is roughly given by 

where the factor of Usn/c accounts for the fact that in the jet 
the photons are free streaming and g ~ 0.15 is a normalisa- 
tion factor. The resulting optical depth to pair production 
is t hen grWer/r, where r i s the radius in the lab frame (see, 
e.g. . iLithwick fc Sarl200lll . Therefore, the lepton density at 
the end of the burst as a result of photon annihilation is 
rigi ~ aTn^r/V, with an associated optical depth along the 
jet of 

1 

T ~ —aTn^±r 




~ 3 X 10"^ , (36) 

for r ~ 10^^ cm. Note that since pair annihilation has been 
ignored, this is only an upper limit. Thus, it may be safely 
concluded that pair production is insignificant during the 
the burst. 

In general, both the jet and the seed photon density and 
temperature will be expected to have radial structure. In the 
jet this is due to competition between the magnetic stresses 
and the Compton drag. In the seed photons this is due to 
the adiabatic cooling of the supernova ejecta. A direct result 
of this structure is that, despite beginning with thermal seed 
photons, the time integrated spectra can have the observed 
broken power law shape where the break energy could indeed 
be interpreted as the temperature of the up-scattered seed 
photons when ry ~ 1, i.e., roughly 2V^ keV ~ few x lO'^ keV. 
For the case where the seed photon temperature has a ra- 
dial power-law dependen ce this has already been explicitly 
shown to be the case bv lGhisellini et al. (2000). However, 
in this scenario, since the density of seed photons depends 
upon distance from the helium star, different orientations 
of the jet with respect to the orbital plane will lead to sub- 
stantial changes in the spectral slopes of the integrated emis- 
sion. Thus, because the compact object is expected to have 
suffered a kick during its birth, the considerable variation 
observed in burst spectral slopes would be expected. 

3.5 Population Statistics 

For a given beaming angle (~ F"^ ) approximate ly lO^^'^F'^ 
bursts occur per galaxy per year JSchmiddbOol . This im- 
plies that the formation rate (TZco-ne) for the progenitor 
systems satisfy 

7^co-Ho > 4 X 10"** galaxy-Vr"' . (37) 



There is a considerable literature which addresses the forma- 
tion rates of neutron star-neutron star binaries and neutron 
star-black hole binaries. Because in both cases it is believed 
that these are produced by the evolution of compact object- 
Helium star binaries, the formation rates of the former place 
lower limits upon the formation rate of the latter. There- 
fore, the progenitor formation rate must be compared to 
TIns-ns — lO^'' to 5 X 10""^ ealaxv~^vr~^ iKaloeera et al' 
^200;^ and 7^NS-BH > 10"" galaxy" Vr"^ iBethe fc Brown 
Il998f) . Hence, there are more than enough presumed com- 
pact object-Helium star products to account for the number 
of 7-ray bursts observed. It should also be noted that the 
formation rate of compact binaries is expected to be con- 
siderably smaller than the formation rate of the progenitor 
binaries due to the possibility of unbinding the progenitor 
as a consequence of the ensuing supernova, and the fact that 
the black hole-black hole binaries have been ignored. 



4 AFTERGLOW 

A consequence of electromagnetic domination in the jet is 
the absence of internal shocks as long as the jet remains 
relativistic. However, as the jet cools, the flow will become 
increasingly hydrodynamic. When F is of order unity, strong 
internal shocks may be expected to develop. At this point 
the internal kinetic energy of the jet can be thermalised, 
producing a moderately relativistic flreball, which can then 
be analysed within the h ighl y successful stan dard fireball 
afterg low model (see, e.g.. lMe szaros 2002; Koni gFfc Grano3 
[200^1 argue that many of the diverse afterglow phenomena 
can be naturally produced if the afterglow is produced inside 
of a pulsar- wind bubble). Due to the velocity gradient at the 
end of the jet, when enough matter accumulates a Compton 
thick head will develop, shutting off further prompt emis- 
sion. This can be expected to occur over the time scale 
for the jet to proceed from the radius at which ry ~ 1 
10^^ cm) to the radius at which F ~ 1 (~ 10^^ cm), which 
for the scenario considered here is approximately 30 — 100 s. 
After this time, the jet will continue to pump energy into the 
growing fireball at its head. Only when the rotational energy 
of the compact object is sufficiently exhausted, or accretion 
ceases, will the jet cease and the fireball expand under its 
own pressure. For a neutron star this also can be expected to 
occur over a comparable time scale as the prompt emission, 
while for a black hole this can continue considerably longer 
(see section I^J. 



5 MODEL IMPLICATIONS 

This model has a number of direct observational implica- 
tions discussed below, many of which have already been de- 
tected. 



5.1 Burst Substructure 

The fact that the ejecta is Thomson thick and inhomoge- 
neous leads to considerable burst substructure. In this case 
seed photons will not be continuously available, but rather 
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will enter the jet in bunches over timescales {5t) associated 
with the inhomogeneity length scales {£), i.e., 



5t 



Is, 



implying 



e ~ 10^ km ~ 10"^i?Hc 



(38) 



(39) 



This will occur when an ejecta clump impacts the jet and 
is subsequently either sheared apart, releasing the entrained 
photons, or forms strong shocks and thus thermalising its 
bulk kinetic energy. 

To lowest order, the rate of supplied seed photons may 
be treated as uniform over St and vanishing otherwise. In 
this CEise, for a single clump, the number of seed photons 
within the jet will evolve according to 



^ - ^iVscod 0<t<5t 
ot ot 



dt 



ot 



St < t 



(40) 



where TV is the total number of seed photons in the clump 
and 



Tb 



FcSt 
2r 



(41) 



which is expected to be ~ 0.3 at the point where most of 
the emission occurs. The solution is trivially found to be 



f T- <t < St 

iV,,ed = AAe""''''' <^ St 

1 St<t 



(42) 



Since nearly every seed photon that enters the jet will be up- 
scattered to 7-ray energies, this also provides the expected 
light curve for a single subpulses 



diV^ _ Tt 



(43) 



where / is the fraction of singly scattered photons and is typ- 
ically of order unity. T his has the fast rise -exponential decay 
structure observed (cf. lNorris et al.lll99(i) . Furthermore, the 
ratio of the exponential decay time to t he (approximately ) 
linear rise time is r^"^ ~ 3, as observed jNorris et alJll99a) . 
As the jet Lorentz factor increases these subpulses would be 
expected to be systematically more symmetric for a given 
inhomogeneity scale in the supernova ejecta (which would 
not be expected to vary considerably). Thus, more energetic 
bursts would be expected to produce more symmetric s ub- 
pulses, which has also been observed (iNorris et alJll993) . 



Many encounters may be approximated by integrating the 
equation 

de £ 



to give 



iVs, 



dA^scat 1 + re/mec2 



-In 



r V e e 



peak 



(45) 



(46) 



and thus e ~ -A~^j in the limit of many scatters. Due to the 
additional path length traversed by these photons, they will 
lag behind the single scattered photons by a time oc Nscat- 
In addition, since the scattered photons are performing a 

1/2 

biased random walk, they will spread in time oc N^^^^. 
gives the following scalings 



This 



and 



Clag ' 



St oc r 



-1/2 -1/2 



(47) 



(48) 



where the former is consistent with the observed anti- 
correlation betwe en the time lags and the overall luminos- 
ity of the burst jBandlll997l) . and the latter is in rough 
agreement with the o bserved relation dlnSt/dlne ~ —0.4 
iFenimore et alJll995^ . The relevant time scales for the time 
lag can be estimated directly by noting that the jet width is 
roughly 



5 X lO-' 



10" cm 



(49) 



and hence the lower energy emission will lag by 
Ascat {r/Tc) ~ 0.2Ascat s. If a majority of the low energy pho- 
tons are produced deeper within the jet (where the Compton 
depth is higher), this time scale can decrease by an order of 
magnitude. In either case this is again consistent with ob- 
servations of burst pulse time lags (.Band.l997i') . 



5.3 Spectral Softening 

As the jet evolves, its axial optical depth will increase and 
the position at which the majority of the emission arises 
from will move outward towards regions of lower Lorentz 
factor. Near the end of the burst the Lorentz factor at the 
7-ray photosphere is expected to be no more than a few. As 
a result, the emission can be expected to soften considerably 
as the burst proceeds. The rate at which this softening oc- 
curs depends upon the radial structure of the jet as discussed 
in section 13.41 and would be expected to vary considerably 
between bursts. 



5.2 Time Lags and Energy Dependent Subpulse 
Widths 

Photons at energies below 2r^T-y will be due to both the 
low energy tail of the seed photon distribution and multiple 
scatters both within the jet and at the jet boundaries. The 
high bulk Lorentz factor of the jet implies that the scattering 
angles will typically be on the order of . Therefore, for 
each encounter 

- —, T ■ (44) 



5.4 Peak— Inferred Isotropic Energy Relation 

lAmati et al.l i2002l) (and more recentlv |Amati2004f) have re- 
ported a correlation between the peak spectral energy and 

1/2 

the inferred isotropic energy of bursts, namely epcak oc EJ^ . 
In the context of the model presented here, this would be 
naturally expected if the dominant factor in the variation 
between bursts was the maximum jet Lorentz factor. This is 
not unexpected if the helium star supernovae are similar, 
implying that scatter in the orbital separation and com- 
pact object parameters, both of which enter most signifi- 
cantly into the determination of F, produces the variability 
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amongst bursts. Then, in terms of the typical seed photon 
energy (esn), the peak 7-ray energy is 

Epeak — 2r^esn • (50) 

From equation 1341 it is clear that 

Liso — Sr** (UT + lis) Usnr^ CX Epcak , (51) 

where the last proportionality holds if the typical ejecta den- 
sities and velocities and the radius at which the jet becomes 
optically thick are similar amongst bursts. This may be triv- 
ially inverted to yield the observed relation. 

5.5 X-ray Flash Characteristics 

For jet Lorentz factors ~ few, the peak of the emission 
would occur in the X-rays, and thus this model provides 
a natural explanation for X-ray flashes. Note that this is 
neith e r a structured nor uniform jet model (cf. iRossi et alJ 
l2002l: iLamb et all l2004l) . In the former, viewmg an az- 
imuthally structured jet from different angles provides the 
peak-inferred isotropic energy relation. In the latter, the 
distribution is in the opening angle of the jet. In both it is 
assumed that there is little scatter in the total energy of the 
burst. In contrast, as seen in equation 13411 with the same as- 
sumptions made in the previous subsection, in the scenario 
presented here the energy released in the form of 7-rays is 
expected to scale as 

Ljot oc P-'' . (52) 

However, in this as well as in the normal bursts, this energy 
is expected to be subdominant relative to the energy in the 
supernova itself. 

Since X-ray flashes share many of the temporal and 
spectral features of 7-ray bursts, including l ying upon th e 
peak-inferred isotropic luminosity relation of lAmatil i2004ri , 
it is tempting to interpret them as simply subluminous 
bursts. Recently, the bolometric energy of the X-ray flash 
XRF 020903 has been measured and was indeed found to 
be similar to typical 7-ray bursts despite h aving a substan- 
tially lower inferred isotropic luminosity dSoderberg et all 

5.6 Late Light Curve Variability 

Despite the large amount of energy released in the super- 
nova, the binary does not necessarily become unbound (in- 
deed, it must not if this is a viable formation scenario for 
double pulsars). Therefore, in systems with a sizable amount 
of ejecta remaining in the vicinity of a nascent compact rem- 
nant, variability due to the occasional accretion by the com- 
panion compact object may be expected after the supernova 
becomes optically thin, and thus late in the light curve. This 
may explain the late time (~ 50 days after the burst) vari- 
ability observed in the residual light curve of SN 2003dh, 
which has a rough time scale of a few days, implying an or- 
bital separation of ~ lO^'^ cm, in agreement with what would 
be expected in this model liviatheson et al...2003»'l . 

5.7 Polarisation 

iLazzati et alJ i2004l) have shown that inverse-Compton scat- 
tering by jets can produce large degrees of linear polarisation 



(~ 100% in some cases). This is a direct result of the rela- 
tivistic aberration of the seed photons in the jet frame and 
is a function of the viewing angle of the jet, nearing unity 
for 6 ^ . Therefore, jet models in which the prompt 
emission is due to inverse-Compton scattering will generally 
exhibit wide variation in the degree of linear polarisation, 
ranging from zero along the jet axis to unity along the jet 
edge. As a consequence, the degree of linear polarisation 
may be used to measure the jet viewing angle, breaking the 
degeneracy between azimuthal jet structure and low lumi- 
nosity in attempts to utilise 7-ray bursts as cosmological 
standard candles. Since, in the model presented here, the 
spectral evolution of each subpulse to lower energy is due 
to multiple scatterings, the polarisation fraction would be 
expected to decrease roughly exponentially with the num- 
ber of scatters and thus the polarisation should be a strong 
function of energy as well. 

Note that as long as the seed photon frequency in the 
jet frame is much greater than the cyclotron frequency, as is 
expected to be the case when the jet is optically thin axially 
(see equation 1221 . the magnetic field will not significantly 
effect the polarisation properties of the prompt emission. 

Presently, the only measurement of the prompt emis- 
sion polarisatio n involved GRB 021206 and was found to be 
80% ± 20% bv lCoburn fc Boggj J2003D (though this claim 
continues to be controversial). However, due to the large 
uncertainty this may be consistent with both synchrotron 
self-Compton models and the model presented here. If this 
is confirmed in future bursts, this would provide strong ev- 
idence for Comptonised jet models. 



6 CONCLUSIONS 

Rough estimates of the consequences of supernovae in com- 
pact object-helium star binaries have been used to suggest 
that such an event is a viable candidate for 7-ray bursts. 
Such events are assumed to occur in the standard double 
pulsar evolution scenario. Within this context a number of 
population calculations have been performed, implying that 
indeed these events are likely to be frequent enough to ex- 
plain the number of bursts seen. This model is capable of 
explaining a number of the burst characteristics, including 
the subpulse light curves and energy dependence, late time 
spectral softening, and the peak-inferred isotropic energy 
relation. In addition. X-ray fiashes and 7-ray bursts are nat- 
urally combined into a single unified theory, as recently sug- 
gested by observations. 

If the binary remains bound, as is necessarily the case 
for the double pulsar formation scenarios, late time variabil- 
ity may be expected. Such variability has been observed in 
the light curve of GRB 030329/SN 2003dh. While no clear 
periodicity is discernible, the typically time scales are con- 
sistent with an orbital separation ~ 10^^ cm as expected by 
this model. Future observations of late time afterglow light 
curves can provide evidence for the existence of binaries in 
7-ray burst progenitors. 

The prediction of significant degrees of polarisation in 
Compton dragged jet models in general, and this model in 
particular, is noteworthy. In the absence of strong struc- 
ture within the jet, this would necessarily lead to significant 
variation in the degree of polarisation, ranging from zero 
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to unity. Predictions for significant degrees of polarisation 
which depend upon viewing angle has significant implica- 
tions for cosmological observations. Even in the absence of 
a detailed understanding of the azimuthal jet structure, po- 
larisation provides a simple way in which to identify the 
viewing angle, removing this degeneracy. Unfortunately, due 
to its considerable uncertainty, the only measurement of 
prompt emission polarisation is as yet unable to differen- 
tiate between synchrotron self-Compton models and Comp- 
tonised jet models. Thus, clearly future polarisation obser- 
vations are required as well. 

Many of the observational predictions of this model 
are generi c to Compton dragg ed jets. Thes e include the 
spect rum (jOhisellini et alJ200(t l. polarisation ijLazzati et alJ 
l200^^j_subpulse lags, and subpulse energy dependence (cf. 
IPado et a l. 2002). To a lesser extent, the Amati relation, 
and thus the X-ray flash characteristics, is also generally ex- 
pected as long as the primary difference among bursts is the 
jet Lorentz factor. 

This model is also very similar to the coUapsar, with 
the obvious difference that in this case the compact ob- 
ject is outside the helium star (which need not be a Wolf- 
Rayet star here). As such, the most immediate feature of the 
coUapsar model, the association with supernova, is present 
here as well. However, it is not possible to simply subsume 
this model into the coUapsar by placing the compact object 
within the helium star; the primary difficulty being that the 
solution to the compactness problem would no longer apply 
due to the considerable supply of MeV photons in the central 
regions of the supernova. As a result the jet would remain 
optically thick until radii two orders of magnitude larger 
than those considered here. Considering the likely fact that 
the supernova ejecta are travelling radially from the com- 
pact object at this point, the mechanism proposed here for 
producing the burst substructure will be unlikely to work. 

Unique to this model is the manner in which the suc- 
cesses of Comptonised jet and coUapsar models are unified. 
Furthermore, it provides a natural class of progenitors, in 
which supernovae are already believed to occur as a result 
of independent observational evidence. Thus, the results pre- 
sented here may be regarded as both a proposal for an new 
7-ray burst mechanism and an investigation into the obser- 
vational consequences of what appears to be an inevitable 
process given the currently observed double neutron star 
systems. 
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